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Summary 


Since the early 1950’s, when Schwartz published his theory of distributions, generalized 
functions have found many applications in various fields of science and engineering. Onte of 
the most useful aspects of this theory in applications is that discontinuous functions can be 
handled as easily as continuous or differentiable functions. This provides a powerful tool in 
formulating and solving many problems of aerodynamics and acoustics. Furthermore, generalized 
function theory elucidates and unifies many ad hoc mathematical approaches used by engineers 
and scientists in these two fields. In this paper, we define generalized functions as continuous 
linear functionals on the space of infinitely differentiable functions with compact support, then 
introduce the concept of generalized differentiation. Generalized differentiation is the most 
important concept in generalized function theory and the applications we present utilize mainly 
this concept. First, some results of classical analysis, such as Leibniz rule of differentiation 
under the integral sign and the divergence theorem, are derived with the generalized function 
theory. The divergence theorem is shown to remain valid for discontinuous vector fields provided 
that all the derivatives are viewed as generalized derivatives. An implication of this is that all 
conservation laws of fluid mechanics are valid, as they stand for discontinuous fields with all 
derivatives treated as generalized derivatives. When the derivatives are written as thjf sum of 
ordinary derivatives and the jump in the field parameters across discontinuities times a delta 
function, the jump conditions can be easily found. For example, the unsteady shock jump 
conditions can be derived from mass and momentum conservation laws. Generalized function 
theory makes this derivation very easy. Other applications of the generalized function theory 
in aerodynamics discussed here are the derivations of general transport theorems for deriving 
governing equations of fluid mechanics, the interpretation of the finite part of divergent integrals, 
the derivation of the Oswatitsch integral equation of transonic flow, and the analysis of velocity 
field discontinuities as sources of vorticity. Applications in aeroacoustics presented here include 
the derivation of the Ivirchhoff formula for moving surfaces, the noise from moving surfaces, and 
shock noise source strength based on the. Ffowcs Williams-Hawkings equation. 

1. Lit ro duct ion 


In the early 1950’s, Schwartz published his theory of distributions that we call generalized 
functions. (See ref. 1.) Earlier, Dirac had introduced the delta function 8(x) by the sifting 
property 

f <30 


/ <p(x)8(x) dx — <p{ 0) 


( 1 . 1 ) 


Dirac recognized that no ordinary function could have the sifting property. Nevertheless, he 
thought of <5(.e) as a useful mathematical object in algebraic manipulations that could be viewed 
as the limit of a sequence of ordinary functions. The Dirac delta function is a generalized function 
in the theory of distributions. Schwartz established rigorously the properties of generalized 
functions. His theory has had an enormous impact on many areas of mathematics, particularly 
on partial differential equations. Generalized function theory has been used in many fields of 
science and engineering. 


To include mathematical objects such as the Dirac delta function into analysis, we must 
somehow extend the concept of a function. The process we use to introduce new objects 
is familiar in mathematics. We extended natural numbers to integers, integers to rationalsj 
and rationals to real numbers. We also extended real numbers to complex numbers. In each 
extension, new objects were introduced in the number system while most properties of the old 
number system were retained. Furthermore, for each extension, we had to think of the new 
number system in a different way from the old system. For example, in going from integers to 



rationals, we view numbers as ordered pairs of integers (a, 6), where 6^0. We identify ordered 
pairs (a, 1) with integer a. The new number system (the rationals) includes the old number 
system (the integers). We must now think of numbers as ordered pairs (a, 6), which we usually 
write as a/b, instead of as a single number a for integers. Similarly, to extend the concept of 
function to include the Dirac delta function, we must think of functions differently. 

We explain in section 2 how to think of functions as functional (i.e., the mapping of a 
suitable function space into scalars). In this way, the Dirac delta function can naturally be 
included in the extended space of functions that we call distributions or generalized functions. 
The usefulness of this theory stems from the powerful operational properties of generalized 
functions. In addition, solutions with discontinuities can be handled easily in the differential 
equation or by using the Green’s function approach. Many ad hoc mathematical methods used 
by engineers and scientists are unified and elucidated by generalized function theory. In fluid 
dynamics, the derivations of transport theorems, conservation laws, and jump conditions are 
facilitated by that theory. Geometric identities for curves, surfaces, and volumes, particularly 
when in motion and deformation, can be derived easily with generalized function theory. In 
section 2 we also define generalized functions as continuous linear functionals on some space 
of test functions. Some operations on generalized functions are defined in this section, as are 
various approaches to introduce generalized functions in mathematics. 

In section 3 we present some definitions and results for generalized functions as well as 
some important results for generalized derivatives, multidimensional delta functions, and the 
finite part of divergent integrals. In section 4 we present various aerodynamic applications 
including derivation of two transport theorems — the interpretation of velocity discontinuity as 
a vortex sheet and the derivation of the Oswatitsch integral equation of transonic flow. The 
aeroacoustic applications include the derivation of the solution of the wave equation with various 
inhomogeneous source terms, the Kirchhoff equation for moving surfaces, the Ffowcs Williams— 
Hawkings equations, and shock noise source strength. All these applications depend on the 
concept of generalized differentiation. Concluding remarks are in section 5 and the references 
follow. 

Many articles and books have been published on the topic of generalized function theory. 
Most of these works have extremely abstract presentations. In particular, multidimensional 
generalized functions, which are most useful in applications, are often treated cursorily in 
applied mathematics and physics books. Of course, some exceptions are available. (See 
refs. 2-7.) Multidimensional generalized functions are relatively easy to learn and use if the 
theory is stripped of some abstraction. To work with multidimensional generalized functions, 
some knowledge of differential geometry and of tensor analysis is required. (See also refs. 8 
and 9.) In this paper, we present the rudiments of generalized function theory for engineers and 
scientists with emphasis on applications in aerodynamics and aeroacoustics. The presentation 
is expository. The intent is to interest readers in the subject and to reveal the power of the 
generalized function theory. Some illustrative mathematical examples are given here to help in 
the understanding of the abstract concepts inherent in generalized functions. 

2. What Are Generalized Functions? 

2.1. Schwartz Functional Approach 

It can be shown from classical Lebesgue integration theory that the Dirac delta function 
cannot be an ordinary function. By an ordinary function we mean a locally Lebesgue integrable 
function (i.e., one that has a finite integral over any bounded region). To include the Dirac 
delta function in mathematics, we must change the way we think of an ordinary function f(x). 
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Conventionally, we think of this function as a table of ordered pairs (x,f(x)). Of course, often 
this table has an uncountably infinite number of ordered pairs. Wf show this table as a curve 
representing the function in a plane. In generalized function theory, we also describe f(x) by a 
table of numbers. These numbers are produced by the relation 

/ OG 

f(x)<f>(x) dx 

-(30 

where the function <f>(x) comes from a given space of functions called the test function space. 
For a fixed function f(x), equation (2.1) is a mapping of the test function space into real or 
complex numbers. Such a mapping is called a functional. We use square brackets to denote 
functional (e.g., F[<f>] and <5 [0] ) . Therefore, a function f(x) is now described by a table of its 
functional values over a given space of test functions. We must first, however, specify the test 
function space. 

The test function space that we use here is the space D of all infinitely differentiable functions 
with bounded support. The support supp <f>(x) of a function <f>(x) is the closure of the set on 
which 4>{x) ^ 0. For an ordinary function f(x), the functional F[<f>] is linear in that, if (f>\ and <p 2 
are in D and if a and (3 are two constants, then 


( 2 . 1 ) 


F[a<p i + l3<t>-2\ - aF[<t>i] + (3F[(f>- 2 ] 


( 2 . 2 ) 


The functional F[<f>] is also continuous in the following sense. Take a sequence of functions {<f> n } 
in D and let this sequence have the following two properties: 

1. There exists a bounded interval I such that for all n, -supp <p n C I ■ 

2. lim <pn\x ) — 0 uniformly for all k — 0, 1, 2, . . .. 
n — kxj 

Such a sequence is said to go to 0 in D and is written <p n 0. Here supp <f> n stands for support 

of <fi n - We then say that the functional F[<f>] is continuous if F[<f > n ] — >■ 0 for <p n — >■ 0. We will 
have more to say in this section about, the space D and why we require the two conditions above 
in the definition of <p n — r 0. 

As an important example of a func tion <f>(x) in D . for a given finite a > 0, we define 


<f>( x;a) 


exp(^) 

0 


( | a; | < a) 
( | .x’ | > a) 


(2.3) 


Note that supp <f>(x a) — [—a, a] and is bounded. We can show that <f>(x',a) is infinitely 
differentiable. Therefore, 4>{x] a) E D. The proof of infinite differentiability at x — ±a is 
somewhat messy and algebraically complicated and we will not belabor this point here. We can 
show that from any continuous function g(x), we can construct another function VK- 1 ') in F> from 
the relation 

4 ! {x) — j g(t)<f>(x — t\a)dt (2.4) 

where the interval [b, c] is finite. The support of xl>(x) is [b — a, c + a], which is bounded. 
The infinite differentiability of i/>(x) follows from infinite differentiability of 4>{x] a). Therefore^ 
tl>(x) E D. There exists an uncountably infinite number of continuous functions. (Consider the 
family of continuous functions sin(aa;), a, f(0,l). This family has an uncountable number of 
members.) It follows from the above argument that there exists an uncountably infinite number 
of functions in space D, so our table constructed from F[<f>] by equation (2.1) representing 
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the ordinary function f(x) has an uncountably infinite number of elements. This fact has an 
important consequence. Two ordinary functions / and g that are not equal in the Lebesgue 
sense (i.e., two functions that are not equal on a set with nonzero measure) generate tables by 
equation (2.1) that differ in some entries. Thus, the space D is so large that the functionals on 
D generated by equation (2.1) can distinguish different ordinary functions. 

We now give an example of a sequence in D such that <f) n 0. Using the. function 

<f>(x',a) in equation (2.3), we define 


<t>n{x) - -<f>{x-,a) 


(2.5) 


This sequence can easily be shown to satisfy the two conditions required for <f> n — 0. We note 
in particular that supp<f> n — [—a, a] for all n. 

Now we define distributions or generalized functions of Schwartz. First, we note that for an 
ordinary function f(x) (i.e., a locally Lebesgue integrable function), the functional F [</>] given by 
equation (2.1) is linear and continuous. The proof of linearity is obvious. The proof of continuity 

requires only that <f> n — r 0 uniformly, which already follows from <p n 0. Remembering that w# : 
are now looking at functions by their table of functional values over the space D and that this 
functional is linear and continuous, we ask if all the continuous linear functionals on spac $ D 
are generated by ordinary functions through the relation given in equation (2.1). We find they 
are not! Some continuous linear functionals on space D are not generated by ordinary functions. 
For example) 

m = m ( 2 . 6 ) 

Proof of linearity is obvious. Continuity follows again from <p n — r 0 uniformly. However, this 
functional has the sifting property that the Dirac delta function requires. As we stated earlier, 
no ordinary function has the sifting property. Therefore, this approach introduces the delta 
function rigorously into mathematics. We define generalized functions as continuous linear 
functionals on space D. The space of generalized functions on D is denoted D 1 . Figure 1 shows 
schematic ally how we extended the space of ordinary func tions to generalized functions. We call 
ordinary functions regular generalized functions, whereas all other generalized functions (such 
as the Dirac delta function) are called singular generalized functions. 

For algebraic manipulations, we retain the notation of ordinary functions for generalized 
functions for convenience.,. We symbolically introduce the notation 6(x) for the Dirac delta 
function by the relation 

5[0]=^(O)= f 8(x)<p(x)dx (2.7) 


Note that the integral on the right side of equation (2.7) does not stand for conventional 
integration of a function. Rather, it stands for 4(0)- We can now use 5(a;) in mathematical 
expressions as if it were an ordinary function. However, we must remember that singular 
generalized functions are not, in general, defined pointwise; they define a functional (i.e., a 
function from our new point of view) when they are multiplied by a test function and appear 
under an integral sign. Thus, when a singular generalized function appears in an expression, it 
is always in an intermediate stage in the solution of a real physical problem. 

More facts about space D in multi dimensions, convergence to 0 in D. and the concept of 
continuity of a functional are appropriate now. The multidimensional test function space D is 
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Space of generalized functions D' 



Singular generalized 
functions 


Real or 

complex numbers 


4*o) m 

Figure 1. Generalized functions are continuous linear functionals on space D of test functions. 


defined as the space of infinitely differentiable functions with bounded support. For example, 
for a > 0, 

‘2 


<f)(x; a) - 




| x | < a ) 

I x I > a ) 


( 2 . 8 ) 


whei 


.. \ 1/2 

22 -c? ) is the Euclidean norm. Other functions in this space can be constructed 
i - 1 ' / 


by using any continuous functiom/(x) and the convolution relation 


V’(x) = / g(t)<p{x - t; a) dt (2.9) 

JQ 

where is a bounded region. The multidimensional generalized function spa®' is defined as 
the space ofcontinuous linear func tionals on the spac <D. In the multidimensional case, a number 
of important singular generalized functions of the delta function type appear in applications. In 
one dimension, the support of6(a;) consists of one point, x — 0. We define the support of a 
generalized function later. In the multidimensional case, in addition fc(x), which has the 
support x = 0, there is also<5(/) with support on the surface /(x) = 0. Section 2.2 contains a 
detailed explanation o©(/). 

We now discuss the definition of continuity of linear functionals on the spa qD. Continuity 
is a topological property. Spac eD is a linear or vector space. It is made into a topological 
vector space by defining the neighborhood of</>(,e) = 0 by a sequence of seminorms. The two 
conditions required above in the definition off) n -3- 0 follow from the conditions used to define 
the neighborhood of — 0. (See refs. 10 and 11.) The definition of continuity of linear 
functionals on space!) can be based on the weak or strong topologies of space!)! (See ref. 11.) 

It so happens that the definitions of continuity based on these topologies are equivalent to the 
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earlier definition of F[<f> n ] — ^ 0 if and only if <j) n -3- 0. Furthermore^ we note that because D is 
a linear space, we can define <fi n ef) if <p n — <f> 5- 0 and because F\<f>\ is linear, we can also say 
that F[<fi\ is continuous if when (f>n—? <f>, we have F \<f>n\ — i- F \<f\. 

We conclude this discussion with one more important fact. Although in this paper we confine 
ourselves to the test function space D , in many applications we should use a different test function 
space. For example, to define Fourier transformation, we should use a test function space S of 
infinitely differentiable functions that go to 0 at infinity faster than |a;| — n for any n > 0 (the 
space of rapidly decreasing functions). Other test function spaces are defined in Carmichael 
and Mitrovic (ref. 10) and in references 12 and 13. Generalized functions on these spaces are 
defined as continuous linear functionals after a suitable definition of convergence to 0 in the 
test, function space is given to get a topological vector space. Note that in all these spaces 
of generalized functions, the important singular generalized functions (such as the Dirac delta 
function) are retained with properties essentially similar to those we study below in space D 1 . 
It can be shown that if A C B and if A and B are two test function spaces used to define 
generalized function spaces A 1 and B 1 , respectively, then we have A C B C B 1 C A 1 (i.e., the 
space of generalized functions A 1 is larger than B 1 ). In particular, D C S C S 1 C D ' , where S 
is thecspace of rapidly decreasing functions defined above. 

2.2. H ow Can Generalized Functions Be Introduced in Mathematics? 

Although Schwartz developed the theory of distributions, like many great ideas in mathe- 
matics and science, the subject has a long history. Synowiec (ref. 14) has stated that evolution 
of the concepts of distribution theory followed a familiar pattern in mathematics of multiple and 
simultaneous discoveries because the appropriate ideas were ‘in the air.’ Several good sources on 
the history of theory of distributions are available. (See refs. 15 and 16.) Therefore, we will not 
present a detailed history here. Also, many different approaches in mathematics can be used to 
introduce and develop systematically generalized function theory. We mention several of them 
here. 

2.2.1. Functional approach to generalized functions. In the functional approach, generalized 
functions are defined as continuous linear functionals. This approach (which we use here) was 
originally introduced by Schwartz (ref. 1) and is the most popular and direct method of studying 
generalized functions. (See refs. 3, 4, and 7.) The operations on ordinary functions such as 
differentiation and Fourier transformation are extended by first writing these operations in the 
language of functionals for ordinary functions, then by using them to define the operations for 
all generalized functions. After the rules of these operations are obtained, the usual notation of 
ordinary functions can be used for all generalized functions. A working knowledge is relatively 
easy to develop with this notation without confusion. Some elementary knowledge of functional 
analysis is needed in this approach. 

2.2.2. Sequential approach to generalized functions. The sequential approach is essentially 
based on the original idea of Dirac in defining a delta function as the limit of a sequence of 
ordinary functions. The approach was originated by Mikusinski from a theorem in distribution 
theory that the space of generalized functions is complete. Therefore, singular generalized 
functions such as the delta function can be defined as the limit of ordinary (i.e., regular) functions 
much like defining irrational numbers as limits of a Cauchy sequence of rational numbers. Many 
good books have been published on this subject. (See refs. 2, 6, and 17.) To define a generalized 
function, the analyst is required to construct and work with a sequence of infinitely differentiable 
functions. Although mathematics only to the level of advanced calculus is involved, the algebraic 
manipulations are technical and laborious. An extension to the multidimensional case also 
appears more difficult than with the functional approach. 
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2.2.3. Bremermann approach to generalized functions. In the Bremermann approach, 
generalized functions of real variables are viewed as the boundary values of analytic functions on 
the real axis. (See refs. 12 and 13.) The Bremermann approach has its basis in earlier works on 
Fourier transformation in the complex plane to define Fourier transforms of polynomials. This 
approach employs some of the powerful results of analytic function theory and is particularly 
useful in Fourier analysis and partial differential equations. A recent book on the subject is by 
Carmichael and Mitrovic. (See ref. 10.) 

2.2.4- Mikusinski approach to generalized functions. The Mikusinski approach is based on 
ideas from abstract algebra. A commutative ring is constructed from functions with support 
on a semi-infinite axis by defining the operations of addition and multiplication as ordinary 
addition and convolution of functions, respectively. This commutative ring has no zero divisors 
by a theorem of Titchmarsh. (See refs. 18 and 20.) Therefore, it can be extended to a field by 
the addition of a multiplicative identity and the multiplicative inverses of all functions. This 
multiplicative identity turns out to be the Dirac delta function. Thfi^ Mikusinski approach gives 
a rigorous explanation of the Heaviside operational calculus and solves other problems such as 
the solution of recursion relations. One limitation of this approach is that the supports of the 
functions are confined to semi-axis or half-space in the multidimensional case. Good sources 
for this approach are Mikusinski (ref. 18), Mikusinski and Boehm# (ref. 19), and an excellent 
expository book by Erdelyi (ref. 20) 

2.2.5. Other approaches to generalized functions. Several other important approaches 
introduce generalized functions in mathematics. One approach is based on the nonstandard 
analysis of Robinson. (See ref. 21.) Nonstandard analysis uses formal logic theory to extend 
the real line by the rigorous inclusion of Leibniz infinitesimals. Many interesting applications 
of this theory, particularly in dynamical systems, are now available* Another more recent 
approach is presented in Colombeau (refs. 22 and 23) and Rosinger (ref. 24). This approach 
uses advanced algebraic and topological concepts to develop a theory of generalized functions in 
which multiplication of arbitrary functions is allowed and it is gaining popularity at present. 
Applications to nonlinear partial differential equations are given by Rosinger (ref. 24) and 
Oberguggenberger (ref. 25). 

3. Some Definitions and Results 

3.1. Introduction 

In this section, some important definitions and results used later are presented. Then, the 
generalized derivative, the multidimensional delta functions, and the finite part of divergent 
integrals are discussed. This paper is application oriented so we are selective about the material 
presented here. We also freely refer to a generalized function by a symbolic or a functional 
notation. 

3.1.1. Multiplication of a generalized function by an infinitely differentiable function. Let 
f(x) be a generalized function in D 1 defined by the functional F[(fi\ and let a(x) be an infinitely 
differentiable function. Then, a(x)f(x) is defined by the rule 

aF[<f>\=F[a<f>\ (3.1) 

Note that aF stands for the functional that defines af. Also, because <f> is in D . so is a<f>. We 
can use this definition to define a(,c)<i(,c). Let <5[<j>] be the Dirac function given by equation (2.6); 
t he n , 

af)[4>] — f)[a4>\ — a(0)#(0) (3.2) 
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Symbolically, this equation is interpreted as 

a(x)6(x) — a(0)6(x) (3.3) 

In the space D 1 , multiplication of two arbitrary generalized functions is not defined; however, 
this statement needs, clarification. Obviously, ordinary functions are also generalized functions 
and any two ordinary functions can be multiplied; thus, they can also be multiplied in the sense 
of distributions. However, multiplication of a regular and a singular generalized function or two 
singular generalized functions may not always be defined. For example, the multiplication of 
6(x) by itself ( i.e. , 5^(a;)) is not defined, neither is f(x)6(x), where f(x) has a jump discontinuity 
or a singularity at x — 0. In applications, experience or inconsistencies in the results occasionally 
show that some multiplications of two generalized functions are not allowed. Sometimes this 
problem can be removed by rewriting the expression such that the troublesome multiplication is 
avoided. For example, difficulties with multiplication of distributions appear if we use the mass 
continuity and momentum equations in nonconservative forms to find shock jump conditions 
(section 4.2). These difficulties can be removed by using the conservation laws in conservative 
forms. To overcome the problem of multiplication of distributions in space D 1 , new spaces of 
generalized functions have been defined. (See refs. 23-27.) Colombeau (ref. 27, chapters f— 3) 
gives an intuitive description of the problem of multiplication of distributions and shows how to 
remedy this problem. 

3.1.2. Shift operator. Let f(x) be an ordinary function and define the shift operator as 
Epf(x) — f(x + h). Then, if F[<f>] is the functional representing f(x) by equation (2.1) and if 
the shifted function Epf(x) is represented by Ef l F\<f>\, we have 

EhF[<f>\ — j /( x + h )<j>( x) dx 

— J f(x)<p(x — h)dx 

- F[E_ h <f>] (3.4) 

This rule can now be used for all generalized functions in D 1 bec ause E_f l <p is in D . For example, 
Epb(x) — 6(x + h ) has the property 

Epb[(f)\ — J 6(x + ft )4>(x) dx 

— 6[ef>(x — ft)] 

= <K~h) (3.5) 

Note that the integral in equation (3.5) is meaningless and stands for the functional Ef,6[f\, 
which in turn is given by 6[E_f l <f>], 

3.1.3. Equality of two generalized functions f(x) and g(x) on an il open set. Two generalized 
functions / and g in D 1 given by functionals F[<p\ and G\<j>\ on D, respectively, are equal on an 
il open set if F[<f>] — G[f\ for all <f> such that supp <f> C 12- For example, d(x) — 0 on open sets 
(0,oc>) and ( — oc>,0). Note that generalized functions are compared only on open intervals. 

3.1. 4- Support of a generalized function. The support of a generalized function f(x) is the 
complement with respect to the real line of the open set on which f(x) — 0. For example, the 
support of 6(x) is the set {0}; that is, the point x — 0. 
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3.1.5. Sequence of generalized functions. A sequence of generalized functions F n [<3] is 
convergent if the sequence of numbers [</>]} is convergent for all ef> in D. For example, let 


8 n (x) — 



0 



UH < 


M > 


(3.6) 


This function is shown in figure 2 and is, of course, an ordinary function. It can be shown that 


lim 8 n (x)—8(x) (3.7) 

n — >oo 

Thus, for the functional <5 n [</>] representing <5 n (,c), 

&n[4>] - j 8 n (x)<p(x) dx (3.8) 

when <f> is in D , we have 

lim 8 n [4>] - (f>(0 ) 
n^oo 

- m ( 3 . 9 ) 


The index in the definition of convergence can be a continuous variable. For example, Fe[<3\ is 

convergent as £ 0 if lim F e \<f\ exists for all <3 in D. 

£ — >0 

The following important theorem characterizes D' and has significant applications. (See 
ref. 7.) 


y 



n n 

Figure 2. Fxampfe of 6 sequence. 
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Theorem: The space D 1 is complete. 

This theorem implies that a convergent sequence of generalized functions in D' always converges 
to a generalized function in D 1 . 

We use this theorem later in this section when we discuss the finite part of divergent integrals. 

3.1.6. Odd and even generalized functions. A generalized function F[<p\ is even if F[<f>(— a;)] 

' ' and odd if F[<f>( — a;)] = — F [</>(,£ )]. For example 6(x) is even and x is odd. 

3.1.7. Derivative of a generalized function. The derivative of a generalized function is the most 
important operation used in this paper. Let f (x) be an ordinary function with a continuous first 
derivative ( i.e. , / is a C 1 function). If /( x) is represented by the functional F[f\ in equation (2.1), 
then we naturally identify its derivative f'{x) with F 1 \<f\ given by the functional 


F'[f\ = 



(3.10) 


Now we integrate by parts and use the fact that <f> has compact support to get 


F'[<p\ = -Jf</>'dx 

= -F[4>'] ( 3 . 11 ) 

Because <f> E D, then / £ f). Thus, Flf'] is a functional on D. We now use equation (3.11) 
to define the derivative of all generalized functions in D 1 . We can keep taking higher order 
derivatives and obtain the following result: 


f (n, M 


( _i)» f L(») 


(n — 1,2,.. .) (3.12) 


We have thus arrived at the following important theorem. 

The orem: Generalized functions have derivatives of all orders. 

We have obtained a very surprising result. Even locally Lebesgue integrable functions that are 
discontinuous are infinitely differentiable as generalized functions. What are the implications of 
this theorem in applications? We address this question about generalized derivatives and their 
applications in section 3.2. First, some examples would be helpful. 

Example 1. The derivative of the delta function F(x) has the property, 

6'[<f>]=-H<f>'] 

— —0'(O) (3.13) 


Symbolically, we can write, 

J 6 l (x)f(x) dx — — <p'(0) 

Note that 6 1 (x) is an odd generalized function. 

Example 2. The H e avis id e function is defined as 


(x > 0) 
(x < 0) 


(3.14) 


(3.15) 


h(x) — 


1 

0 
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or in functional notation, 


H[(f>]— / <p(x)dx 

J 0 


(3.16) 


This function is discontinuous at x — 0. To define the generalized derivative, we use 

equation (3.11) as follows: 


H '[</>] -H[</>' ] 



= m 

= S[<f>] (3.17) 


Symbolically, we write 

h'(x) — 6(x) (3.18) 

Note the use of the bar over h 1 to signify generalized differentiation because h 1 (x) — 0 where 
now h 1 stands for the ordinary derivative. 

We give one more important characterization of space D 1 (ref. 7) known as the structure 
theorem of distribution theory. 

The orem: Generalized functions in D' are generalized derivatives of a finite order of continuous 
functions. 


For example, we note that the Dirac delta function is the second generalized derivative of the 
continuous function 


/(*) = 



(x > 0 ) 
(x < 0 ) 


(3.19) 


3.1.8. Fourier transforms of generalized functions. We now work with the space of rapidly 
decreasing test functions S . (See sec. 2 .1, the last paragraph.) In this space the Fourier transform 
of each test function is again in S. (See refs. 2, 4, 6, and 7.) We define the Fourier transform of 
an ordinary function il>(x) as 

/* oc> 

$(£)= / f(x)e 27rix t dx ( 3.20) 


Let f(x) be an ordinary function that has the Fourier transform /(£) (e.g., let / be square 
integrable on (— oo, oo)). Then for il>(x) in S, the Parseval relation is 



f(x)f(x) dx 



f(x)f( a;) dx 


(3.21) 


If now F[V’] is identified with f{x), then we should identify F[V’] with /(^). However, 
equation ( 3.21 ) is actually the relation 


F[rf>]=F$] (3.22) 

We use this relation as the definition of the Fourier transform of generalized functions in space 
S ' . For example, 

^ ^ COG 

<S[t/>] — 5[^i] — ^(0) — / dx (3.23) 

J ^oc 
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The last integral is the functional generated by the function 1 so that 

H) = 1 (3.24) 

Thus, the Fourier transform of the Dirac delta function is the constant function 1. 

We will not discuss this subject further because we do not use Fourier transforms extensively 
in this paper. We note, however, that if ip i s in D , then ip is not necessarily in D and 
equation (3.22) is meaningless in D 1 . Therefore, we must change the test, function space from D 
to S. Another method of fixing this problem is to use the Fourier transforms of functions 
in space D as a new test, function space D. The Fourier transformations of functions in 
D' are now continuous linear functionals on space D. These generalized functions are called 
ultradistributions. (See ref. 28.) 

3.1.9. Exchange of limit processes. One of the most, powerful results in generalized function 
theory is that the limit processes can be exchanged. For example, all the following exchanges 
are permissible: 



(3 .25a) 
(3.25b) 
(3.25c) 






1 

y 






-271 —71 

0 

71 

-1 

2k 







Figure 3. Periodic function with jump discontinuity of 2(— 1 )" at x = me, n = 0, ±1, ±2, . . •• 
Also, by equation (3.25b), we have 


dj_ 

dx 


i <3C> . 

d 4 

dx ^ (2m + 1 ) 7T 

777 =0 V ' 


oc 

E 

777=0 

(30 


dx (2m + 1 )7r 


— cos(2?7i + 1) 

• ' '■ * TT 


777 =0 


sin(2?77-|- l)a; 


sin(2?77 + 1 )a; 


x 


From equations (3.28) and (3.29), we conclude that 

2 (30) (30 

— COs(2?77 + 1 ).e = ( — 1 ) n 6(x — 77 7r) 

m~ 0 7i—— oo 


(3.29) 


(3.30) 


The series on the left is divergent in the classical sense. Nevertheless, such a result is often useful 
in signal analysis. Another important application of exchange of limit processes is in obtaining 
the finite part, of a divergent integral. (See sec. 3.4.) 

3.1.10. Integration of generalized functions. We say that G[(f\ is an integral ofF[<(>] if 


G'[<p] - F[ef>] 

(3.31) 

For example, we can easily show that the Heaviside function is 
function because 

an integral of the Dirac de It 

-HW] 


= <K«) 


,-e-, 

II 

(3.32) 

Let K[f\ be a generalized function such that 


K 1 [<j)\ = 0 

(3.33) 
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for all <f> £ D. Then, if G[(f>\ is an integral of F[<fi], it follows that (G+ K)[(f>] — G[(f>] + K[<f>\ is 
also an integral of F[(f\. References 7 and 29 show that the only solution of equation (3.33) in 
D ' is 

A-M =/«,(,)«, (3.34) 

where c is an arbitrary constant (i.e., K[f\ is a constant distribution). This result corresponds 
to the classical indefinite integration of a function 

J f(x)dx — g(x) + c (3.35) 

We use the same notation symbolically for all generalized functions. For example, we write 

JsM d , = hM +c ,3.36) 

where h(x) is the Heaviside function. Note that the integral on the left of equation (3.36) is 
meaningless in terms of the classical integration theories. 

3.2. Generalized Derivative 


The generalized differentiation concept is quite important in generalized function theory; 
this section focuses on it and gives some very useful results for applications. Indeed, the results 
themselves, rather than the mathematical rigor used in deriving them, are of interest in this 
paper. As before, a bar over the differentiation symbol denotes generalized derivatives if there 
is an ambiguity in interpretation. For example, we use df/dx, f'( x ), 3f/3x and d f /3x; dx j 
to denote generalized derivatives of ordinary functions, but we do not use a bar over 8 1 (x) and 
38(f) [3xj because it is obvious that these derivatives can only be generalized derivatives because 
8(x) and 8(f) are singular generalized functions. 

3.2.1. Functions with discontinuities in one dimension. Let f(x) be a piecewise smooth 
function with one discontinuity at with a jump at this point defined by the relation 


A/ = f(x o+) - f(x o_) 


(3.37) 


We want to find the generalized derivative of f(x). Let cf> be in D and let ,cq be in the 
support of 4>(x). Then if F\<f>\ is the functional representing f(x) by equation (2.1), we have for 
supp 4> — [a, 6], the result 


F'm = —F u'l 


/ (x)3\x) dx 


fX 0 - fh 

/ + / f(x)<f) l (x)dx 

f a ^^0+ 


f'(x)<f)(x)dx+ [f(XQ + ) - f(XQ-)]<f)(XQ) 


f l (x)<p(x) dx + Af <P(xq) 


(3.38) 
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We have performed an integration by parts to get to the last step. We have also used the fact that 
4>{a) — 4>{b) — 0 in the integration by parts. Noting that <f> Gq) — 6 \<f> (x + ^q)] = E_ A . 0 <5[/(,c)], 
where E— X q is the shift operator, we write equation (3.38) symbolically as 

./Vi = f'(x) + AfS(x-x 0 ) (3.39) 

One question is the use of f f (x) compared with the ordinary derivative f f (x). Let us study 
equation (3.38). The functional F'[<f\ corresponding to f'(x) indeed has retained the memory of 
the jump A / on the right side of the equation. Symbolically, f 1 (x) can be integrated over [c , it?] , 
where c < xg < x, to give the result 

/ X 

f'(x) dx + /(c) + A fh(x - x 0 ) (3.40) 

Thus, we have recovered the original discontinuous function. We note, however, that 

/(*) / jy , (x)dx + f(c) (3.41) 

because the memory of the jump A / is not retained in f'(x) but is retained in f'(x ). If a 
function f (x) has n discontinuities at x;, i — 1 — n with the jump Af; at x; defined by 


Af;- fix;. i (3.42) 


f'(x) = fix) + £ A f; 6(X - x • i (3.43) 

/= 1 

This equation is the first indication that when we work with discontinuous functions in 
applications, the proper setting for the problem is in the space of generalized functions. In 
particular, if an integral method, such as the approach that uses the Green’s function, is used 
to find the solution, essentially no significant changes to algebraic manipulations are needed in 
finding discontinuous solutions provided we stay in the space of generalized functions. Again, 
we will have more to say about this later. (See sec. 3.2.3.) 


3.2.2. Functions with discontinuities in multidimensions. Let us now consider the function 
/(x), which is discontinuous across the surface g(x) — 0. Let us define the jump Af across 
g — 0 by the relation 

A/ =f{g = 0+) -f{g = 0-) (3.44) 

Note that eg — 0+ is on the side of the surface g — 0 into which Vg points. We would like to 
find df/dx;. To do this we use the results from section 3.2.1 as follows. Let us put a surface 
coordinate system ( V v?) on g — 0 and extend the coordinates to the space in the vicinity of 
this surface along normals. Let v? — g be the third coordinate variable that is well defined by 
the function g in the vicinity of this same surface. We note that / in variables u and u is 
continuous, but it is discontinuous in variable iA Therefore, we have 


W_ = d£_ 

Out Out 


(4=1,2) 


■df_ 

duf 


df_ 

du 3 


+ Af6(J) 


(3 .45a) 
(3.45b) 
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In equation (3.45b), we used equation (3.39). Thus, using the summation convention on index 
j, we get 


Of Of duJ 

Ox 1 duJ Ox; 


df OuJ 
dui dx-i 


+ A / 



9 / 

-d^ + Af d^ 6{9) 


We can write this in vector notation as 


V/ = V/ + AfVg6(g) 


(3.46) 


(3.47) 


In section 3.2.3, we discuss how to interpret 6(g) when g — 0 is a surface. We can similarly 
define generalized divergence and curl as follows: 

V f = V f + Vg • Af 6(g) (3.48a) 

Vxf = Vxf + VjxAf«(y) (3.48b) 

The rigorous derivation of both these results requires some knowledge of the invariant definition 
of divergence and curl in general curvilinear coordinate systems. (See refs. 8 and 9.) We can 
combine the above three results by using * for the three operations such that 

V *f = V*f + Vg * Af 6(g) (3.49) 


3.2.3. Ordinary differential equations and Green's function. We give a few simple results 
here. One important question discussed in connection with integrals of generalized functions is 
the solution of 

./'m-i = 0 (3,50) 

in D 1 . It can be shown easily that the only solution of this equation is the classical one (refs. 7 
and 29) 

f(x) = C (3,51) 

where C is a constant. However, the solution of the equation 

xf(x) - 0 (3,52) 

which is not a differential equation, is 

f(x) - C 6 (x) (3,53) 

To get this solution, some simple results from the generalized Fourier transform are used. (See 
ref. 29.) Taking the Fourier transform of both sides of equation (3.52), we get 

= 0 (3,54) 

df 
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Therefore, after integration of equation (3.54), we have 


/(£) = C 


(3,55) 


By taking the inverse Fourier transform of both sides of equation (3.55), we get equation (3.53). 
From this result, the solution of 

xf'(x) - 1 (3,56) 

is found as 

f(x) — In | | + Cq + C 2 h ( x ) (3.57) 

where C\ and C 2 are constants and h(x) is the Heaviside function. The solution C 2 h\x) comes 
from the fact that the generalized function fjx) satisfying the equation 

xf'(x) - 0 (3,58) 


is, from equation (3.53), 


f'(x) - G 2 6 (x) 


(3,59) 


Thus, the solution of the homogeneous equation (3.58) is the integral of this function 


f(x) - C\ + C 2 h(x) 


(3.60) 


Let us now consider a second order linear ordinary differential equation with two linear and 
homogeneous boundary conditions (BC) as follows: 

&u. — A(x)u" + B(x)u' + C(x) — f(x) 1 (a- 6 [0, 1]) 

BCi[u] = 0 > (3.61) 

BC 2 [«] = 0 J 

Let us also assume that we know u is a function and u n is Lebesgue integrable so that u n — u n 

and u — v! . Suppose a function g(x,y) exists, the Green’s function, such that 

«(•£)= / f(y)g{z,y)dy (3.62) 

Jo 

Because u G C'\ then (u — (u by continuity of u and u . We know we can take ( into the 

integral in equation (3.62) but not ( because g(x,y) may not belong to CJ . Therefore, using ( x 

to indicate that derivatives in i are with respect to x, we get 

I u — (u 


= ?x f(y)g(x,y)dy 

Jo 

- / f(y)?xg(x,y)dy 

Jo 

= /(*) 

from equation (3.61). We see that i x g{x,y) must have the sifting property 

'.rgix.y) - 6(x - y) 


(3.63) 


(3.64) 
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Because the BC’s are linear, we also have 


BC[u] = [ f(y)BC x [g{x,y)]dy (3.65) 

JO 

Therefore, other conditions on g(x,y) are 

B C i A . [sr (it? , y)] - 0 (3.66a) 

BC2a;[fif(a;, y)] - 0 (3.66b) 

where the x in the subscripts of the BC’s indicates that g(x,y) in the variable x satisfies the two 
boundary conditions. 

From equation (3.64) we conclude that, because £ x is a second order ordinary differential 
equation, g(x,y) must be continuous at x — y and dg/dx must have a jump discontinuity at 
x — y. The reason is that if g(x, y) has a discontinuity at x — y, the first generalized derivative 
with respect to x will give a b(x — y) by equation (3.39). A second generalized derivative would 
give 6'(x — y) in the result. But because 6 1 (x — y) is missing on the right of equation (3.64), 
g(x,y) cannot be discontinuous at x — y. Assuming that g(x,y) is defined by 


v{x,y) 


9^,y) 

9 2 {?, y) 


(0 < x < y) 
(y < x < 1) 


equation (3.64) means that 


1 //I - ' 1 - 0 

9 x {y,y) - 9. 2 (y,y ) 

dg _ dg, 1 

-7 I (y,y) - -r L (y,y) - 

ox ox A(y) 


(3.67) 


(3 -68a) 
(3.68b) 

(3.68c) 


Note that equation (3.68a) is the same as £ x g — 0 used above. This equation means that g ^ and 
g in variable x are solutions of the homogeneous equation £u — 0. Equation (3.68b) expresses 
continuity of g at x — y and equation (3.68c) gives the jump of < dg/dx at x — y. To get 
equation (3.68c), we note that 


?x9 - ?x9 + A {y) 


dg ‘2( \ dg U \ 

-M-(y> y) - -M-(y> y) 

ox ox 


6(x - y) 


= Mv) 


dg r 
dx 


{y,y) 


dg 1 

dx 


{y,y) 


6(x - y) 


- - y) 


(3.69) 


The last delta function follows from equation (3.64). Equation (3.68c) follows from the fact that 
the coefficient of 6(x — y) in the expression after the second equality sign must be equal to 1. 
The Green’s function is now determined from equations (3.66) and (3.68a— c). 


3.2.4. Leibniz rule, of differentiation under the integral sign. We want, to find the result, of 
taking the derivative with respect, to variable a in the following expression in which A, B, and 
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/ are continuous functions and B( a) > .4(a) for a E [a, 6]. Thus 


d 

E( a) — — / f(x,a)dx 
da .1 


Let us define the function H(x, a) as follows: 

H(x,a) — h[x — A(a)]h[B(a) — r-] 


(3.70) 


(3.71) 


where h(x) is the Heaviside function. The function H(x, a) — 1 when A(a) < x < B( a) and 
H( x,a) — 0 otherwise. Using H(x, a), we can write E(a) as 


d f 00 

E(a) — — / H (r 1 , a)f(x, a) dx 

da J x 


r3C> /9H f a. u df '\ w 
T—f + H — ) dx 


da ' 


da 


(3.72) 


We have 


(x, a) — —A'(a)h[B(a) — a;] 5 [a; — .4(a)] 

da 

+ B l (a)h[x — .4(a)]<i[i?(a) — a;] 

= — A' ( a)5[a; — .4(a)] + B 1 (a)6[B(a) — a;] (3.73) 


Note that we have used 

h[B(a) — a;] 5 [a; — .4(a)] = h[B(a) — .4(a)]<i[a: — .4(a)] 

= $[a;-A(a)] (3.74) 

because B(a) — .4(a) > 0; thus, the Heaviside function is 1. Similarly, we do the same as in 
equation (3.74) for the second product of the Heaviside and the delta functions in equation (3.73). 
Using equation (3.73) in equation (3.72) and integrating with respect to x, we get the Leibniz 
rule of differentiation under the integral sign, 

fB(n) q* 

E(a) — I — — (x , a) dx B (a) f[B(a) , a] — A(a) f[A(a) , a] (3.75) 

JA(a) Oa 


3.3. Multidimensional Delta Functions 


In multi dimensions, 5(x) has a simple interpretation given by 

J d(x)6(x) dx — < 56 ( 0 ) 


(3.76) 


Thus, 


<5(x ) = <S(a;i )<S(a: 2 ) . . .6( x n ) 


(3.77) 


where x = (a:p x%, . . . ,x n ). In this section, we confine ourselves to three-dimensional space. Of 
interest in applications are 8(f) and 8'(f) where / = 0 is a surface in three-dimensional space. 
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We can always assume that / is defined so that | V/| = 1 at every point on / = 0. If / does not 
have this property, then _/q — //|V/| does. Thus, redefine the surface. 


3.3.1. Interpretation of 6(f). Consider the integral 

I - j dx 


(3.78) 


Assume that we define a curvilinear coordinate system (t/^, t/ 2 ) on the surface f — 0 and extend 
these variables locally to the space near this surface along local normals. Let iP — f which, 
because |V/| = df / dip — 1, t/ 3 is the local distance from the surface. Thus, / = iP — constant 
0 is a surface parallel to / = 0. Of course, we assume iP is small. From differential geometry 
(refs. 8 and 9), we have 

dx = ^ j ( w'* , u 2 , iP ) dip du 2 dip (3.79) 

where g^(iP,iP, iP) is the determinant of the first fundamental form of the surface / = iP — 
Constant. Using equation (3.79) in equation (3.78) and integrating with respect to iP gives 


I = 


1 

< 1 

2 

3\ 

X 1 

V ’ 

« , 

u ) 

xl 

(w 1 , 

P 2 , 

°I 


6 ^t/ 3 ^ \j C J( 2)(u'*, iP,iP) dv} dip- dip 


.1 n ,.‘2 


'/= 0 


<p(x)dS 


That is, I is the surface integral of <f> over the surface / = 0. 


(3.80) 


3.3.2. Interpretation of P{f). We want to interpret 


I - J <P(x)fi'(f) d* 


1 2 3 

x | u , u , u 


P ( u 3 ) \ I g (2 ) ( u 1 , u 2 , u 3 ) du 1 du 2 du 3 


(3.81) 


Here we have used the coordinate system (t/1, t/ 2 , tr 3 ) defined above. Integrating the above 
equation with respect to iP gives 


I = 


diP 


^( x )\/d( 2 )( u1 U /2 W 3 ) 


i / :1 =0 


dip dip 


(3.82) 


Again, from differential geometry, we have 


3 

d^\/d { -2) 


I u‘ ,iP ,iP) — —2 Hj{^iP,iP,iP'J \Jyp2) ( u 


'. u -’.- p ' 


(3.83) 
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where H j- stands for the local mean curvature of the surface / = — Constant. Taking the 

derivative of the integrand of equation (3.82) and using the result of (3.83), we obtain 


I = - 


d<(> 

du? 


x ^ u 1 , u 2 , 0 ^ \j 9{ 2)( u '* ! « 2 , 0) d u 1 d i? 


+ J '2H j , u 2 , 0^</> x^t/lj t/ 2 ,!)^ \j 9{ 2 )( u '*, « 2 , 0) dv} du 2 


'/= 0 L 




dS 


(3.84) 


where d<f>/dn is the usual normal derivative of 0 . Intuitively, the appearance of the term 2 Hjffi 
in the integrand is not at all obvious. This appearance is a clear indication of the importance 
of differential geometry in multidimensional generalized function theory. 


3.3.3. A simple trick. We have already shown that <fi(z)8(z) — f(0)8(x). By taking the 
derivative of both sides of this equation, we get 


<p , (z)8(z) + <p(z)8 , (z) — <fi( 0)8' (z) 


(3.85) 


Obviously, the right side is simpler than the left side. Let us consider the expression 


E = <t>{x)6{f) = f 




(3.86) 


where again we have used the coordinate system 
know that 


1 2 3 

I xl u , u , u 


(t/1, t/ 2 , tr^) defined in section 3.3.1 above. We 


<t> 


|x| u } , t/ 2 ,0 



(3.87) 


We use the notation <p(x) for <fi [x (u^ , t/ 2 , 0) ] ; that is, <p(x) is the restriction of <p( x) to the support 
of the delta function that is the surface / = 0. We note that d<j)/dn — (9/9u^)^[x (t/l,t/ 2 ,0)] 
= 0. Using <f>(x), we can write E in two forms: 


E — <p(x)8(f) (First form) 

E — <fi(x)8(f) (Second form) 

Is there an advantage of using the second form compared with the first form? The answer is 
yes! Let us take the gradient of E for the two forms in equation (3.88). Thus, 

VE = Vf 6(f) + <f>{x) V/ 6\f) (First, form) 

VE — V-jf 8(f) + <fi(x) V/ 8'(f) (Second form) 

Here, V -jf i s the surface gradient of </>(x) on / = 0. From equation (3.84), we note that 
in the integration of 8'(f) in the first form, the term d<f>/dn cancels a similar term in the in- 
tegration of V0 8(f). In the second form, because df/dn — 0, d<f>/dn does not appear in 
the integration of 8'(f) and obviously is also absent in the integration of V2 4>8(f). Therefore, 



(3.88) 
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algebraic manipulations are reduced. It is, thus, expedient to restrict functions multiplying 
the Dirac delta function to the support of the delta function. Note carefully that functions 
multiplying 8 1 (x) cannot be restricted to the support of 8\x); that is, <p(x)8 l ( x) ^ <p(0)8 l (x). 


3.3.4- The divergence theorem revisited. Let it be a finite volume in space and let </>(x) be a 
C 1 vector field. Let us define the discontinuous vector field <f> i(x) as 


4>(x) (x e 9) 

0 (x ^ 9) 


(3.90) 


Let the surface f — 0 denote the boundary 99 of region 9 in such a way that n = V/ points to 
the outside of 99 and |V/| = 1 on / = 0. We have 


= V •</>!+ A^-n 6(f) 

= V-0i-^(x)-n 6(f) (3.91) 

We note that A </>i = </>i(/ = 0+) — 4>\{f — 0 — ) = — 0). Integrating over the unbounded 

three dimensional space, we get 


dx\ 


dx\ dx -2 dx 3 




dx% — 0 


(3.92) 


Similarly, we get zero for integrals of 9<jq 2 / an d 9<jq 3/9^3, where 9] is the ith component 
of (j)\ . Therefore, 

j V • <f>i dx = 0 (3.93) 

Now, the integration of the right side of equation (3.91) using equation (3.80) gives 


V ■ 4> dx — <f> n dS — 0 

Jq Joq 


Here we have used the fact that, from equation (3.90), 


V • = 


V • <9 
0 


(x e 9) 
(x 0 9) 


(3.94) 


(3.95) 


Also, we define <p n — <f> • n. Equation (3.94) is the divergence theorem. 

We note that equation (3.93) is valid if <p\ has a discontinuity across the surface k — 0 within 
9 as shown in figure 4. Equation ( 3.94) is therefore valid if V- <f> in the volume integral is replaced 
by V • <f>, where the only jump of <f> in the generalized divergence comes from the discontinuity 
on k — 0. That is, we write 

V- <9= V •</> + A<f>-n'6{k) (3.96) 

where n/ = Vfc is the unit normal to k — 0. Equation (3.94) can now be written 



V ■ 4>dx — 



(3.97) 
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Vk = n' 



Figure 4. Control volume Q intersecting surface of discontinuity of vector field <f> used for deriving generalLed 
di verge nc e th eorem . 

which, by using equation (3. 96), we can also write as 

[ V • 4 > dx — [ <f> n dS — [ A <f> n idS (3.98) 

Jn JdQ Js k 

where A <ft n i — A <f> ■ n 1 and Sj, is the part of the surface k — 0 enclosed in region (See fig. 4.) 

The divergence theorem is used in deriving conservation laws in fluid mechanics and physics 
in differential form. The fact that it remains valid for discontinuous vector fields, as shown 
in equation (3.97), implies that such conservation laws are valid when all the derivatives 
are interpreted as generated derivatives. Thus, the jump conditions across the surface of 
discontinuities are inherent in these conservation laws as shown in section 3.4. This interpretation 
of conservation laws eliminates the need for the pillbox analysis of jump conditions. 

3.3.5. Product of two delta functions. We have said earlier that the product of two arbitrary 
generalized functions generally may not be defined. Here we give the interpretation of the 
product of two multidimensional generalized functions for which multiplication is possible. Let 
f — Q and g — 0 be two surfaces intersecting along a curve T as shown in figure 5. Assume 
V/ = n and Vg — n/, where |n| = |n/| = 1. We want to interpret 


/ = J <p(x)6(f)6(g) dx 


(3.99) 
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On the local plane normal to the T-curve, defineu = /, u — g. and u — 1. where f is the 
distance along the T-curve. Extend t/ 1 and u 2 to the space in the vicinity of the plane along a 
local normal to the plane. We have 


dv} du 2 du^ 1 



(3.100) 
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Vg = n' 



Figure 5. Integration of 6(f) 6(g) for two intersecting surfaces / = 0 and g = 0. 


where sin# = |n x n , |. Using equation (3.100) in equation (3.99) and integrating the resulting 
integral with respect to w' and we get 

I = f LA du 1 di? du 3 = / — IdT (3 .101 ) 

J sm# V / V / Jr q sm# 

9=0 

This result is useful in applications. (See sec. 4.3.) 

3.4. Finite Part of Divergent Integrals 

The finite part of divergent integrals is important in aerodynamics. The classical procedure 
for finding the finite part of divergent integrals appears ad hoc and leads to questions about the 
validity of the procedure. First, could the appearance of divergent integrals in applications be 
the result of errors in modeling the physics of the problem? Second, will the method lead to a 
unique analytical expression or do different analytical expressions lead to equivalent numerical 
results? The generalized function theory clearly answers these questions. 

Let us first examine the function /(a;) = In | .c | , which is locally integrable. The ordinary 
derivative of this func tion is 

——hi I a; I = — (3.102) 

dx x 

which is not locally integrable over any interval that includesa; = 0. We know, however, that 
as a generalized function, ln|a;| has generalized derivatives of all orders. What is the relation of 
the generalized derivative of In | .:c | to the ordinary derivative/^ a;) = 1/a?? 

Let us work with F[<f> ] representing ln| a? | as follows: 


m = 


In |a? | </>( a?) dx 


{<P e D) 


We have, using the definition of generalized derivative, 



ln|a:|/ / (a:) dx 


(3.103) 


(3.104) 
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Figure 6. Function h. E (x) n-ivl in defining finite part of divergent integral J[^>(x)/x]dx. c > 0; rv > 0. 


We need some integration by parts to get the term ]/ x in the integrand of equation (3.104). 
However, this integration cannot be performed because x is not locally integrable. We solve 
this problem by using a new functional depending oir, the limit of which isF 7 ^] as follows. 
Let h,s(x) be a function defined below for some constants > 0 and a parameter e > 0. Thus, 


h c {x) — 


(s<x<as) 
(Otherwise ) 


(3.105) 


This function is shown in figureb. Then it is obvious that 1 i^.c | can be written as the limit of 
an indexed generalized function as follows 


li in //.ff.r i In x 

£-— n 


ln|a 


(3.106) 


Note that if we define Fl[<f)\ as 


Fl[<p\ ——I /ig- ( a? ) In | it? | ( a; ) dx 


(3.107) 


then we have from the completeness theorem of D' (sec. 3.1.5) 


limFiM =F'[ 


(3.108) 


The function he (a? )ln|it? | has two jump discontinuities at a; = —e and x — as. We can either 
apply the classical integration by parts to equation (3.107) by breaking the real line into two 
intervals or by using the generalized derivative 


F' £ [<f>] — J — [/)■£■ (a;) ln|;t? |]0( a? ) dx 


(3.109) 


Here we are integrating oveisupp (j) and we do not worry about the: terms coming from the limits 
of the integral in the integration by parts because^ = 0 at the limit points. 

We now take the derivative of the term in square brackets in equation (3.109): 


d 

— [h £ (x)ln\x\\ — — lns6(x + s) + ln(o5gr)5(a? — as) 

dx x 


(3.110) 
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Here and below, we have used the result that <f>(x)6 (x — ^'0 ) — 4> ( ) 6 (x — £(]). Thus, after 
using equation (3.110) in equation (3.109) and integrating with respect to x, we have 


f hc(x) 

Fl[<p\ — — in e <f>( — e ) + ln( ae )<fi(as ) + j <p{ x) dx 

J x 

f hAx\ 

= <^(0)lna + / <f>(x)dx + o(e) 


(3.111) 


where o(e) stands for terms of order e and higher. Now from equation (3.108) we have 


F '[</>] = 


f hc-ix) 

— <t>( 0)ln a + lim / <f>(x) dx 

e— 0 J x 


— <f)( 0)ln a + lim 



<{>{z) 

X 


dx + 




x 


(3.112) 


We can show that the limit of the integral on the right of equation (3.112) exists. If now a —l, 
then In a- = 0 and 


F'[4>] - lim 


<Ma) 


dx + 


<Kx) 


dx 


(3.113) 


which is known as the Cauchy principal value (PV) of the integral. But a — 1 need not be taken 
and equation (3.112) is numerically the same as equation (3.113). The above limit procedure is 
called taking the finite part of a divergent integral. 


What have we achieved? Over any open interval that does not include x — 0 we have 


—hi la; I = — (3 .114) 

dx x 


but when x — 0 belongs to the open interval, then the classically divergent integral must 
be interpreted such that the functional F'\<f>\ corresponding to (d/dx) In |.t; | is recovered. As 
the above simple function demonstrates, more than one different analytical expression for the 
procedure can bf used to find the finite part of a divergent integral. However, all the expressions 
are numerically equivalent. We define the principal value of l/x as 


PV 



d 

— ——In la; 

dx 


(3.115) 


Thus, when x — 0 is in the interval of integration of 1/a;, the finite part of the divergent integral 
must by taken to get the numerical value of F'\<f>\, where F[<f>] is given by equation (3.103). Note 
that the term regularizing a divergent integral is also used in mathematics. The procedure given 
here corresponds to the canonical regularization of Gel ’f and and Shilov. (See ref. 7.) 

What is the use of this procedure in applications? Suppose we have reduced the solution of 
a problem to the evaluation of the expression 


u 




<P(y) ln k -y\dy 


(3.116) 
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where x E -h Let us assume that we know that the integral is continuous as a function of x so 
that d/dx can be replaced by d/dx and taken inside the integral. We get 


f d 

u(x) - / <p(y)—hi\x - y\dy 
JQ dx 


= / 0(y )P V 

J Q 


1 


x - y 


dy 


(3.117) 


which is interpreted as the finite part of the divergent integral by the procedure defined earlier. 
We remind the readers that the procedure will result in exactly what equation (3.116) would 
give had we been able to perform the integration analytically. Also, assuming that <f> — 0 at the 
boundaries of 0, an integration by parts of the first integral in equation (3.117) would give 

u(x) — — (p 1 (y)hi\x — y\dy (3.118) 

which is also a legitimate result if this integral exists. The problem is that often in applications, 
equation (3.116) is an integral equation for the unknown function <f>(x), which has integrable 
singularities at the boundaries of the interval f2. Thus, the above integration by parts is invalid 
and, in any case, the integral equation (3.118) is divergent. Therefore, the only choice left is 
the integral equation with the principal value of l/(x— y), which is a well-known kernel in the 
theory of singular integral equations. 

We now give an advanced example in three dimensions with a surprising implication in the 
numerical solution of an integral equation of transonic flow which we will discuss in section 4. 
Let us consider the integral 

« x) =*[m i y (3.119) 

dx f Jn r 

r' 2 - (aq - yi) 2 * + {x‘ 2 - y 2 f + (£3 - ys ) 2 (3 .120) 

where Q is a region in space and x 6 O. In this problem </>(x) is a function and is the 
unknown of the aerodynamic problem. Assuming that the integral is a function in x, we can 
replace d z /3x^ with 3 /dx^ and take the derivatives inside the integral 

Hx)= L m id^) dy 

=L«’'-k(j ) d7 13 1211 

We use generalized differentiation rather than ordinary differentiation because the latter will 
result in a divergent integral. Note that 

3 _ ?*i 

dy 1 \r J 

3 2 /l\ 3 r 2 — r 2 

dy 2 V7 ?* 5 



(3.122a) 
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where r\ — a;i — y\ . Because r\/r^ is integrable, we write 


w_ 

dy\ 


M A 

(H) 

1 dyi 

VrV 


(3.123) 


and we proceed to find the finite part of the divergent integral in equation (3.121). 

Let /( y,x, £") = (?*i, ?*2 , r. 3 ) — e — 0 be a piecewise smooth surface enclosing the point 

y = x where ?q; = Xj — y;, i — 1—3 and g is a homogeneous function of order 1; that is, 

g ( ar\ , ar-j , aiq) = ag(i q,?* 2 , 7 * 3 ). This condition assures that the surface g [ar\ , 0 - 7*2 j ar 3 ) — £" — 0 
corresponds to g ( ?q , 7*2 , 7 * 3 ) — e/a — 0 for a ^ 0. Thus, all the surfaces g — e — 0 correspond 

to various values of a that are similar in shape. From the homogeneity of g, it follows that 

/( y,x, 0 ) — g (?*i, ?* 2 , 7 * 3 ) = 0 consists of a single-point y = x. For example, for a sphere with a 
center at y = x and radius e, we have 


/( y , x, £■) = \J r\ + ?| + ?| - f = 0 


(3.124) 


In addition, we assume 'Vyf — n, where n is the local unit outward normal to the surface. Let 
/ > 0 outside and / < 0 inside this surface, respectively. We introduce the function hs{ y) by 
the relation 


My) = 



(/ > 0 ) 
(/ < 0 ) 


(3.125) 


Now, we define the required generalized derivative in equation (3.123) by the relation 


dy'i 


\ • cT 

= lim — — 


/ e— 0 dy\ 



lim 


M (/) + - 


-My) 


(3.126) 


where n\ is the component of n along the yj-axis. Therefore, I(x) can be written 

/(x) = lim / 1 -^-<l)(y)dS 

, 3 2_ 2 

+ lim / ,r l r h £ (y)<p(y)dy (3.127) 

e — 0 Jq 7*° 

where we have used equation (3.80) to integrate S(f) in equation (3.126). 

LI sing a Taylor series expansion of </>(y) y = x, we find that 


lim 



Mp-</>(y) dS - aj<f>(x) 


(3.128) 


where a j- is a constant depending on the shape of the surface / = 0. For example, for the sphere 
given by equation (3.124), we have 


47T 



(3.129) 
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If we take the surface / = 0 to be a circular cylinder with its axis parallel to the yi-axis such 
that the base radius is s and its height is 7 , £- ^> 1 , then 

a f = 4tt (3.130) 


Equation (3.127) is thus written 


/(x) = aj-^(x) + lirn f 3? 1 ? h € {y)<f>{y)dy (3.131) 

Numerically, / (x ) is the same regardless of the shape of / =0. Because (3rj — r^) /r^neary = x 
takes both positive and negative values, the shape of / = 0 as e —? 0 affects the value of the 
integral in the summation process. This effect is similar to a well-known result for conditionally 
convergent series, which can be made to converge to any value by rearranging the terms of the 
series. The term a in equation (3.131) compensates for the change in the value of the 
volume integral when / = 0 is changed so that /(x) is numerically the same. 

What is the implication of the above result in applications? In practice, the volume 
integration is performed numerically. The volume integral has a hole enclosing y = x whose 
boundary surface is given by / = 0. The value of aj must, therefore, correspond to the grid 
system used in the volume integration. If the hole is rectangular, which is often the case, then 
neither of the above two a j-'s in equations (3.129) and (3.130) is appropriate for the problem. 

One question remains unanswered. When does the appearance of a divergent integral imply 
anything other than the breakdown of the physical modeling? The answer is when we have 
wrongly taken an ordinary derivative inside an improper integral. Such a step can make the 
integral divergent and is caused by thfiwrong mathematics (improper procedure) rather than the 
wrong physics. Thus, the analyst should always check the cause of the appearance of divergent 
integrals in applications. Because in classical aerodynamics, the inappropriate mathematics 
generally causes the appearance of divergent integrals, the finite part of divergent integrals must 
be used. 

4. Applications 

4.1. Introduction 

In this section we give some applications in aerodynamics and aeroacoustics that show 
the power and the beauty of generalized function theory. We use the results of the previous 
sections here. Many areas of aerodynamics and aeroacoustics can use generalized function 
theory, especially because the approach is almost always more direct and simpler than other 
methods. In addition, for many problems involving partial differential equations, no alternate 
method is available for finding a solution. Below is a partial list of applications of generalized 
function theory in aerodynamics, fluid mechanics, and aeroacoustics: 

Aerodynamics and fluid mechanics 
Derivation of transport theorems 

Derivation of governing conservation laws (such as two-phase flows) 

Derivation of jump conditions across flow discontinuities, velocity discontinuity as a vortex 
sheet 

Derivation of the governing equation for boundary element or field panel methods 
Subsonic, transonic, and supersonic aerodynamic theory 
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Ae roacoustics 


Sound from moving singularities 

Derivation of the governing equation for the boundary element method 

Derivation of the Kirchhoff formula for moving surfaces 

Study of noise from moving surfaces using the acoustic analogy 

Identification of new noise generation mechanisms and their source strength (such as shock 
noise ) 

In addition, in both aerodynamics and aero acoustics, generalized function theory can help in 
the derivation of geometric identities involving curves, surfaces, and volumes, particularly under 
deformation and in motion. 


4.2. Aerodynamic Applications 


We give here four applications that have been previously derived by other classical methods. 
The method based on generalized function theory, as expected, is much shorter and more elegant. 
Other examples in aerodynamics are presented by De Jager. (See ref. 5.) 


4-2.1. Two transport theorems. We give two results here that are used in the derivation of 
conservation laws. We want to take the time derivative inside the integral 


d f 

1 - — / Q(x,t)dx 

at jQ(t) 


(4.1) 


where it(t ) is a time-dependent region of space and Q(x,t) is a C'd function. Let us assume the 
boundary dit(t ) of it is piecewise smooth and is given by the surface f — 0 such that / > 0 in 
it. Assume also that V/ = n/ where n/ is the unit inward normal to the surface. Suppose we 
can ascertain that the integral in equation (4.1) is continuous in time. Then, we can replace 
d/dt with d/dt and bring the derivative inside the integral. We write 


I - J t ' J Hf)Q{^,t)drx 


^-6(f)Q(x,t) + h(f)^ 


dx 


df 


3Q 


'dQ(t) 


d t Q{ X’t )dS + J m m 


(4.2) 


where h(f) is the Heaviside function. Here we have used equation (3.80) to integrate 6(f) in the 
second step above. We can show that 


dj_ 

dt 


- - V = v n 


(4.3) 


where v n i and v n are the local normal velocities in the direction of inward and outward normals, 
respectively. Thus, 

I - [ v n Q(x,t) dS + [ ^-dex (4.4) 

JdQ(t) JQ(t) dt 

This equation is the generalization of the Leibniz rule of differentiation of integrals in one 
dimension. 
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For the second result, we want to take the time derivative inside the following integral by 
assuming again that the integral is continuous in time and that Q is a C'^ function. Thus, 

* = -£/ Q(x,t)dS (4.5) 

We first convert the surface integral into a volume integral 

1 = ^ j H f)Q{*J) d * (4.6) 

Here / = 0 describes dit(t ) and V/ = n, where n is the unit outward normal. Also, note that 
Qis the restriction of Q to / =0 as explained in section 3.3. Therefore, 


I = 


If 9Q ' 


dx 


(4.7) 


We now must use the results of section 3.3 to integrate 6'(f) and 6(f). However, df / dt — —v n 
and this function is restricted already to / = 0. Thus, 


d 

dn 


i>n,Q(x , t) 


= 0 


(4.8) 


Using equations (3.80) and (3.84), we get 


I = 



~dQ 

dt 


2 v n HfQ(x,t) 


dS' 


(4.9) 


where H j- is the local mean curvature of dit(t) . What is dQ / dtl We have the following result, 
assuming that Q is no ni impulsive, 


__ dQ dQ_ 

dt dt n dn 


(4.10) 


Derivation of equation (4.9) by other methods is not trivial. 


4-2.2. Unsteady shock jump conditions. These conditions are usually obtained by the pillbox 
analysis. Wf present a method here based on generalized function theory. We have said that 
the conservation laws such as the mass continuity and the momentum equations are valid as 
they stand if we replace all ordinary derivatives with generalized derivatives. We derive here the 
jump conditions from these two conservation laws. Let k(x,t) — 0 describe an unsteady shock 
surface. Let Vfc = n, where n is the unit normal pointing in the downstream direction. We 
denote this downstream region as region 2 and the upstream region as region 1. We define the 
jump A Q in any parameter by 

AQ — [Q]% — [Q\i (4.11) 
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Wing 


Wake 


Figure 7. Thin wing in incompressible, irrotational flow with wake. 


where the subscripts 1 and 2 refer to the upstream and downstream regions, respectively. 
Applying the rules of generalized differentiation to the mass continuity equation, we have 


f L + y > u) = f7 +y ■()>») 


+ 


(5 k 

Ap dt + A(/m) ' 11 


6(k) = 0 


(4.12) 


where p is the density andu is the fluid velocity. The sum of the first two terms on the right of 
the first equality sign is the ordinary mass continuity equation and is 0. The coefficient ^k) 
must also be 0. Thus, 

-Apv n + A (pu n ) - A \p (u n - v,j)] = 0 (4.13) 

where v n — —dk/dt is the local shock normal velocity and« n = u • n is the local fluid normal 
velocity. This expression is the first shock jump condition. 

Tile momentum equation in tensor notation using the summation convention gives 


1 A (wii) + A_ = !(„„,) + A ( „„,„. ) + A 

Qt l> 7 ),.. 1 J’ fl„- rtf l> fl„- 1 J> fl„- 


OK r 


dx; dt 


dxj ■ Ox; 


+ 


3 k 


A (pit;) — + A [pujUj) n j + A pn; 


6(k) = 0 


(4.14) 


where p is the pressure. The sum of the three terms after the first equality sign is the ordinary 
momentum equation. The coefficient of 6(k) must be zero; therefore, 


A [pu;(un - ('?>)] + A pn; — 0 (4.15) 

This expression is the second shock jump condition. We can derive a similar result from the 
energy equation. 

Note that had we used the mass continuity and momentum equations in nonconservative, 
form, we would have been faced with ambiguities of multiplication of generalized functions. 
This problem is discussed in detail by Colombeau. (See ref. 27.) In that reference, the remedy 
for the removal of these ambiguities is discussed from intuitive and mathematically rigorous 
aspects. 


4.2.3. Velocity discontinuity as a vortex sheet. Let us consider a thin lifting wing in forward 
flight in an incompressible fluid as shown in figure 7. It can be shown that the velocity field can 
be idealized as irrotational (i.e.y x u = 0) where uis the fluid velocity. However, a velocity 
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Figure 8. Diagram used in deriving Oswatitsch integral equation of transonic flow. 


discontinuity can occur on the wing and on the wake. In this cas57 x u / 0 and the velocity 
discontinuity over the wing and the wake gives the vorticity distribution 

V x u = V xu+nx Au<5 ( k ) 

= u x Au.8(k) 

= r 8{k) (4.16) 

where k(x,t) — 0 describes the wing and wake surfaces andVfe = n, the local unit normal to 
these surfaces. Here we define the vorticity distributionT = n x Au. Note that again we define 
Au = [u ]-2 — [u]i, where n points into region 2. The Biot-Savart law gives the velocity field, 

u(x)= m^8(k)dy= [ ^dS (4.17) 

J i Jk=0 r 


where 


r 


x y 


r 


4.2.4. An integral equation of transonic flo w. To derive the Oswatitsch integral equation of 
transonic flow (refs. 30 and 31), consider a thin wing with shock waves in transonic flow moving 
with uniform speed along the aq- axis as shown in figure. 8. Let u be the perturbation velocity 
along the aq-axis. The governing equation for this flow parameter in nondimensional form is 

9 1 

v 2 7hf = 0 <“*> 

For simplicity, we assume that the airfoil, the shock surfaces, and the wake surface are all 
specified by k(x) — 0. We set up this problem in generalized function space by converting 
the derivatives in equation (4.18) to generalized derivatives. We again define a jump iaa or 
ir by A(-) = [-(2 — [ •] 1 , where n = Vfe points into region 2. For the airfoil itself, we define 
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[«] i = u 2 . = 0 bec ause the airfoil is a closed surface. Thus, 


-2,„ 2 


— o 1 W u 

V u - 


2 d i. 


= A 


' du 

d 

( “ 2 11 

3n 

dn\ 

(tJJ 


+ V • < A 


-ni 


6(k) 

6(k) 


(419) 


where n = (np n% n%) is the unit normal to the surface k — 0 and n | = (np 0,0). Note that 
on the right side of equation (4.19) we dropped the sum of two terms, which by equation (4.18) 
is 0. 

To get an integral equation, we use the Green’s function of the Laplace equation, which is 
— 1/47 rr, and treat — < d 2 u 2 /dx 2 as a source term to obtain 


47ru(x) = — 


i d ' 


2 dx'{ 


1 9 

-u (y) dy 


1 

du 

d 

( u ' 2 \ 

■0 r 

_ On 

dn\ 

V 2 ). 


dS 


V-/ — fun — -7—111 ) dS 

Jk = 0 r V 2 


(4.20) 


Now if we; bring the derivatives inside the first volume integral, which is over the unbounded 
space, we must use the finite part of the divergent integral introduced in section 3.4, equa- 
tion (3.119). Taking Q( y) = « 2 (y) in equation (3.119), from equation (3.131) we have 


3, 


— / -u 2 (y)dy- aju 2 (x) + lini 


3 ?* 2 


5 h £ (y)u‘ 2 ( y) dy 


(4.21) 


The last integral in equation (4.20) is 

V-/ - fun - ill ) dS = 

Jk = 0 r V 2 


>k = 0 ?* 


— ( u cos 9 — — n\ cos 9\ ) dS (4.22) 


where cos^ = r • 11, cos 9\ — 1/ n\ r • lip and r = (xj — yi )/?*. Our job is finished. The integral 
on the right side is convergent. Substitute equations (4.21) and (4.22) in equation (4.20). The 
result is the Oswatitsch integral equation of transonic flow. Further approximation is possible, 
but we stop at this point. This derivation is much shorter and more direct than the original one. 
(See refs. 30 and 31.) 

4 . 3 . Aeroacoustic Applications 

In this section, we give four examples for the linear wave equation. Even for this equation, the 
use of generalized function theory leads to important and useful results. Before the examples, 
we give some standard forms of the inhomogeneous source terms appearing in aeroacoustic 
problems. These follow: 

□ 2 $ = Q(x,t) 

□ 2 $ = Q(x, *)$(/) 
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(4.23a) 

(4.23b) 



(4.23c) 


n 2 * = 2-[Q{x,m.f)] 

at 

□ 2 $ = V- [Q(x,t)6(f)\ (4.23d) 

□ 2 $ = Q(x, *)/»(/)$'(/) (4.23e) 

n 2 $ = Q{x,t)6{f)6(f) (4.23f) 


In these equations, f(x,t) — 0 is a moving surface, usually assumed a closed surface. An 
open surface, such as a panel on a rotor blade, is described by / = 0 and f(x,t) > 0, where 
f(x,t) — f(x,t) — 0 describes the edge of the open surface. (See fig. 9.) Also, we denote the 

Heaviside function as . In equation (4.23e), note that Q is the restric tion of Q to f — 0. 

The solutions of the above Equations have been given in many publications of the author and 
coworkers. (See refs. 33-36.) We give only a brief summary here. 


The Green’s function of the wave equation is 


G( y, r; x,t) 



( T < *) 

(r>t) 


where 


r 

9 - r-t + - 

c 


(4.24) 


(4.25) 


In this equation, (x,t) and (y,r) are t h < observer and the field (source) space-time variables; 
respectively. The speed of sound is denoted by c and r — |x— y|. The two forms of the solution 



Fgure 9. Defmitbn of open surface by relations / = 0, / > 0. Edge is defined by / 
geodesic normal. 


/ = 0 and v is the unit inward 
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of equation (4.23a) are 


(4.26) 


47r$(x,/) = J^[Q]ret,dy 

and 

4tt ®(x,t) = f - —f Q(y,T)<m (4.27) 

J —oo t T JQ(t) 

where the subscript ret stands for retarded time t — r/c. The surface f2(r) is the sphere 
r — c(t — t) (i.e., the sphere with center at the observer x and radius c(t — r) with the element 
of the surface denoted by di}) . The two forms of the solution of equation (4.23a) are known as 
the retarded time and the collapsing sphere forms of the solution, respectively. 

The solution of equation (4.23b) can also be written in several forms. (See refs. 33 and 34.) 
We give two forms here. For a rigid surface /( x, t ) = 0, let M r — M r be the local M ach number 
in the radiation direction. Then 

4tt<F(xT) = / i dS (4.28) 

7/ =0 [r|l-M r |J ret 

Note that to get this, equation, the formal Green’s function solution, which is 

47t$(x,/)= [ -Q(y,r)6(f)6(g)dy dr (4.29) 

Jr 

is integrated as follows. First introduce a Lagrangian variable // on and near the surface f — 0 
such that the Jacobian of the transformation is unity. Note that we have; y = y (ij,t) and 

r - |x - y(i], r)| (4.30) 

Next let t — r g , which gives dg/dr— 1 — M r . Integrate equation (4.29) next with respect to g 
and finally integrate 6(f) by the method of section 3.3 to get equation (4.28). 

A more interesting method of integrating the delta functions in equation (4.29) is to let r — r g 
and integrate with respect to g . The integration gives 

47T$(x, /) = j ~[Q( y, T )] ret 5( F) dy (4.31 ) 


where F(y;x,t) - [/(y,r)] ret = f[y,t - r/c]. Note, however, that even if |V/| = 1 by 
definition, we have |VT| 1 in equation (4.31). We will, therefore, give a slight modification of 
equation (3.80). In the following integral, assume |V/| 1, then 

/= f <f>(x)6(f)dx= f (4.32) 

J Jf= 0 |V/| 

This result applies to equation (4.31). It is easily shown by differentiation that 


1 

| VF| = (l + M'l - 2 M n cos dj 7 = A 


(4.33) 


where M n — v n / c , v n — —df/dt is the local normal velocity on / = 0, and cos 9 — n r is the 
cosine of the angle between the local normal to / = 0 and the radiation direction r = (x — y )/r. 
Using equations (4.32) and (4.33) in equation (4.31), we get 

dS (4.34) 

ret 
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47T$(x,/ ) 


1 Q{ y,r) 
r A 




Figure 10. Influence surface for observer on propeller surface, forward Mach number = 0.334, helical Mach 

number = 0.880. (See ref. 44.) 

where c?E is the element of the surface area at F — 0. Note that for supersonic surfaces, the 
condition M r — 1 produces a singularity in equation (4.28). The us* of equation (4.34) removes 
this singularity in most cases. 

To visualize the surface E: F — 0, let the surface / = 0 move in space. Construct the inter- 
section of the collapsing sphere?* = c(t — r) for a fixed (x, /) with the surface / = 0. The surface 
in space that is the locus of these curves of intersection is the E-surface or thfqifluence surface 
for (x, t ). Given (x,t ), this surface is uniquebecause the sphere r — c(t — r) has center at x and 
?* = 0 at t — t. Given r <t, because (x,t) is fixed, the sphere is specified and/(y, r) = 0 is also 
specified. Therefore, the intersecting curve, if it exists, is specified. Figure 10 (ref. 32) shows the 
E-surface for a rotating propeller blade. This figure indicates that the E-surface is dependent 
on the motion and the geometry of the surface / = 0. A singularity in equation (4.34) may exist 
when A = 0; however, we will not address the singularity problem here. Such a problem can 
occur for supersonic propeller blades with blunt leading edges. That situation should be avoided 
because of excessive drag problems. 

The solutions of equations (4.23c) and (4.23d) can be related to that of (4.23b). For example, 
the solution of equation (4.23c) is 


dt J f =Q [?* |1 - M r 


Now d/dt can be brought inside the integral using the relation 


d _ 1 d 

dt 1 — M r dr 


(4.35) 


(4.36) 


Note, however, that r — |x — y(? 7 ,r)| so that dr/dr ^ 0. (See refs. 33 and 34.) Similar 
manipulations can be performed for the form of solution of equation (4.23c) based on equation 
(4.34), but it is better to work with the source terms of equation (4.23c) before using the Green’s 
function approach. (See section 4.3.3.) 
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The solution of equation (4.23e) is by far the most difficult of the problems considered here. 
We first simplify the algebraic manipulations by defining / such that V/ = v where v is the 
unit outward geodesic normal to the edge. The geodesic normal is tangent to the surface / = 0, 
/ > 0 and is orthogonal to edge f — f — 0. (See fig. 9.) The formal solution of equation (4.23e) 
is 


4tt<4>(xT) = f ^-Q(y,T)h(f)6 l (f)6(g) dydr 
.! r 

- J r )]ret h(F)6'(F) dy 


(4.37) 


where, as before, F — [/] re t and we define F(y; x, /) = [/(y,r)] re t- Let N be the unit normal to 
the surface F — 0. We can show that 


N = 


n — M n v 

A 


(4.38) 


Equation (4.37) is now of the form of equation (3.81). Again, because |VT| = A ^ 1, we must 
give the modification to equation (3.81) here. In this case, we have for |V/| ^ 1 

Next, using F in place of / here, we get from equation (4.39) 


IV/ 1 


2 H f 4>) 

+ ww\ db 


(4.39) 


47r<f>(x, t) 



/[Qjret \ 2 H F [Q\ ldt h(F)] 

A^l^— ^ ) J + ^ 


(4.40) 


where Hp is the local mean curvature of the E-surface given by F — 0. Note that 

JLh{F) = K.VFb(F) (4.41) 

so that we must integrate this delta function in equation (4.40). Using a curvilinear coordinate 
system on the E-surface, we can show that 



<ft(x) |V F|5( F) dS 



F= 0 


(4.42) 


where 9 1 is the angle between N and N = VF/|VF|. Also, dL is the element of length of the 
edge of the E-surface given by F — F — 0. The final result of manipulations of equation (4.40) 
based on equation (4.42) is 


47r < f>(x,f ) 



F> 0 


1 d 

A W 



2H F [Q]ret) 


+ 


rA 2 




dS 


[Q] ret col 


rA 2 


-dL 


(4.43) 


F = 0 
F= 0 
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Note that we have defined N as the unit normal to F — 0 and we have 


N = 


■ - M„ r 


(4.44a) 


A 

A=|VF| (4.44b) 

cos 9' — N ■ N (4.44c) 

Because Q is the restriction of Q to f — 0, we have dQ/dn — 0. In equation (4.43), we find the 
normal derivative of [Q]ret &st: 


d 

In 


[Q]ret — 


~dQ- 

i 

+ - 

' dQ 

dN 

— N • r 

c 

dr 

_ 

ret 

- 


(4.45) 


ret 


Using equation (4.38), we get 


N = — [(1 — M n cos 9 ) n — M n sin 9t J 


(4.46) 


where t] is the unit vector along i lie projection of r on the local tangent plane to / = 0. 
Therefore, after using dQ/dn — 0 we get 


dQ 


M n sin# 


dQ 


dN a dti {4Al) 

where d/dt\ is the directional derivative of Qalong tj. In this case, we no longer need restriction 
of Q to / = 0 because dQ/ dt i — dQ/dt For dQ/ dr in equation (4.45), we must use a relation 

similar to that in equation (4.10). The curve F — F — 0 in equation (4.43) is generated by the 
intersection of the collapsing sphere g — 0 and the edge curve / = / = 0 of the open surface 
/-»•/.> »• 

We next consider equation (4.23f). The formal solution is 

47 t $( x ,/) = I ^Q(y,T)6(f)6(f)6(g)dydr 




(4.48) 


This equation is similar to equation (3.99) except that |VF| ^ 1 and |VF| ^ 1. In this 
case, sin# in equation (3.100) is replaced by |VF x VF|, which by definition is Aq. Therefore, 
equation (4.48) gives 

47r$(x,/ ) — [ — [ ® 

Jf = n r 


>F = 0 

f = n 


Ao 


d.L 


(4.49) 


J ret 


We now give four applications. 


Jf.3.1. Low-son's formula for a dipole in motion. A dipole is an idealization of a point force. 
A point force in motion is described by the wave equation 


□V = {/,;(* )% - s (*)]} 

ox i 


(4,50) 
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where p 1 is the acoustic pressure, f; is the component of the point force, and s(t ) is the position 
of the force at time t . The formal solution of equation (4.50) is 


4tt p'(x,t) - J ^p~6 [ y - s(r)]%) dydr 

Let us integrate the above integral with respect to y. We get 


(4,51) 


47rp / (x,i 


(4,52) 

where 

r* — x — s( r) 

(4,53a) 

and 

* 



r 

g - r-t H 

c 

(4,53b) 

Now let t — r g* and note that 

-r— — 1 — M r 

OT 

(4,54) 


where M r = s t/c is the Mach number of the point force in the radiation direction. Integrate 
the resulting equation with respect to g* to get 


4:irp l (x,t ) 


d 

dx-i 


fii T ) 

r 1 1 — M r | * 


(4,55) 


where r* is the emission time. The solution of g* =0 has only one root if the point force is in 
subsonic motion. The derivative in equation (4.55) can now be taken inside the square brackets. 
Thfbresulting equation is a formula given by Lowson (ref. 37) that is useful in noise prediction 
of rotating blades where the dipole sources can by assumed compact. 


4-3.2. Kirchhoff form ula for moving surfaces. In the 1930s, Morgans published a paper in 
which he derived the Kirchhoff formula for moving surfaces. (See ref. 38.) The derivation of this 
formula was based on classical analysis and was lengthy. In 1988, Farassat and Myers gave a 
modern derivation of this result based on generalized function theory. (See refs. 39 and 40.) The 
derivation is short and avoids the use of four-dimensional Green’s identity and the associated 
difficulties of dealing with surfaces and volumes in four dimensions. We present the basic idea 
behind this modern derivation here and refer the readers to reference 39. 

Assume that the surface in motion on which conditions on <f(x,t) are specified is given by 
f(x,t). This surface can be deformable. Assume that <f> satisfies the wave equation in the 
exterior of the surface / = 0, which is the region defined by / >0. Now extend <f> to the entire 
unbounded space as follows: 


<p{x, t) 


<p(x,t) 

0 


(/> 0 ) 
(/< 0 ) 


(4,56) 


Clearly, <f> satisfies the wave equation in the unbounded space. However, <f> has discontinuities 
across / = 0 that appear as source terms of the wave equation. Note that the jumps in 
<fi and its derivatives depend on corresponding values for <f> because A<f> — <p(f — 0+) and 

a df/dt ~ dfj dt (/ = o+ ). 
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Applying the rules of generalized differentiation to <ft , we get 


□ > = 


+ l -M n ^8{f) - I|- [M n <f>6(f)\ - V • [<f>n6{f)] 


(4.57) 


where M n — v n /c and v n — —df/dt is the local normal velocity of the surface / = 0. As before, 
we have assumed V/ = n, the local outward unit normal to / = 0. The three types of source 
terms on the right of equation (4.57) are of the standard types given in equations (4.23a— f ). The 
solution for a deformable surface is given by 


4irip(x,t ) = 


D(S) L r ( 1 — M r)\ T * 


dv} du 2 + / 

J D(S) 


' ^ 2 ^ 2 )- 
r (1 - M r ) _ , 


1 c\ 

du du 


(4.58) 


where D(S) is a time-independent region in u^u^-space onto which the surface / = 0 is mapped. 
The determinant of the coefficient of the first fundamental form is denoted In this equation 

r* is the emission time of the point (t/I ,v?) on the surface / = 0. The expression E\ depends on 
<ft, <p n , V‘ 2 <f> (surface gradient of <f>), and the kinematic and geometric parameters of the surface 
/ = 0. The expression E-j depends only on the kinematic and geometric parameters of the 
surface / = 0. (See ref. 39.) 


4-3. 3. Noise from moving surfaces. Let an impenetrable surface / = 0 be in motion such 
that / > 0 outside the body and V/ = n, the unit outward normal. Let us assume that the 
fluid is extended inside this surface with conditions of undisturbed medium (i.e., density po and 
speed of sound c). We know that the mass continuity and momentum equations are valid when 
the derivatives are written as generalized derivatives. Let us extract only the contribution of 
discontinuities across f — 0 and leave the effect of all other discontinuities (such as those across 
shock waves) in these equations. The mass continuity equation gives 

Op — 

+ V • (pu) = - [p - pa)v n 6( f) + puj( f) 

- PQVnH f) (4.59) 


where v n — —df/dt is the local normal velocity of / = 0 and we have used the impenetrability 
condition on this surfac^. which is u n — v n . The momentum equation gives 

“7 /"'A 1 + — ( /'"I"./ + P ij) = P /i (4.60) 

where Pjj = E,q + (p — po)6jj is the compressive stress tensor and E,q is the viscous stress 
tensor. Now we take the generalized derivative of both sides of equation (4.59) and d/dxj of 
both sides of equation (4.60), subtract the latter from the former, and finally subtract d p/dx'j 
from both sides to get 

□V = [%*</)] - £-[ P «’>y(/>] + f[«l»„«(/)] (4.61) 

where p' — c^(p — pq). Here Tj-j is the Lighthill stress tensor. Now we have added h(f), the 
Heaviside function, on the right side to indicate that Tjj 0 outside the surface / = 0. This is 
the Ffowcs Williams-Hawkings (FW-H) equation. (Sefqref. 41.) Note that in the far field, p' is 
the acoustic pressure. 
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The source terms of equation (4.61) are of the standard types in equations (4.23a— f). For 
a surface in subsonic motion, the solution for surface sources involving the Doppler factor is 
most appropriate for numerical work. (See refs. 33 and 34.) For supersonic surfaces such as 
an advanced propeller blade on which \M n \ < 1 everywhere, a different solution based on the 
E-surface must be used. We show here briefly how this can be done. In applications, we need to 
calculate the sound from an open surface such as a panel on a blade. We, therefore., define such 
an open surface by / = 0,/ > 0 with the edge defined by / = / = 0 as before. The assumptions 
concerning the gradients of / and / at the beginning of this section hold here. We are interested 
in the solution of equations of the types 


□ 

□ 




d 


p - - 


Ox; 


POVnM'mf) 

pnih('f)Hf) 


(4.62a) 

(4.62b) 


where h(f) is the Heaviside function. Note that we have approximated P ; j in equation (4.61) by 
p6; j where p is the surface (gauge) pressure. To derive solutions for equations (4.62a) and (4.62b) 
suitable for supersonic panel motion, we write 


d 

dt 


PQVnh(f)H.f) 


dv n 

po w H'm.f) 


- PQVnHf)^(f) 


- Pa v n v,,Hf)6(f) 


(4.63) 


where v v is the velocity of the edge in the direction of the geodesic normal. A similar operation 
can be performed on the right of equation (4.62b) such that 


d 

dx; 


pnih(f)S(f) 


= -V 


pnh(f)6(f) 


= -ph{f)6\f) + 2pH f h{f)6{f) 


(4.64) 


where H j- is the local mean curvature of the surface / = 0. The source terms of the right of 
equations (4.63) and (4.64) are of the types in equations (4.23a— f). (See refs. 35 and 36.) 


4.°3.4. Identification of shock noise source strength. The first term on the right of the 
FW-H equation (4.61) is known as the quadrupole source. As mentioned in the derivations, 
the discontinuities in the region / > 0 (i.e., outside the body) contribute source terms after 
generalized differentiation is performed. If a shock wave described by the equation k(x,t) — 0 
exists on a rotating blade, then the quadrupole term gives surface sources on the shock, the 
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strengths of which are determined as follows. Let us take the generalized second derivative of 

T ij h y 


dTjj 

dxj 

9 2 t- • 
u >./ 

dx i dx j 


3T; 


Ox.: 


■ + AT ijn' i 6(k) 


T; ; 


>./ 


Ox: dx , 


dT- • 


dx , 


(4.64a) 

(4.65b) 


where n! — Vfe is the unit normal to the shock surface pointing to the downstream region. The 
last two terms on the right of this equation are shock surface terms that, are of monopole and 
dipole types, respectively. The first, term on the right, of equation (4.65b) is a volume term that 
is familiar in Light. hill’s jet. noise theory. In the rotating blade noise problem, this term primarily 
reflects nonlinearities other than turbulence. Farassat and Tadghighi (ref. 42) conjectured that 
the shock surface terms contributed relatively more than the volume term in equation (4.65b). 
Preliminary calculations have supported this conjecture. (See ref. 43.) 

The interesting aspect, in the above result, is that the shock source strength is obtained 
purely by mathematics. Without, the use of the operational properties of generalized functions, 
the identification of shocks as sources of sound and the determination of the source strength 
would be rather difficult.. Other mechanisms of noise generation can also be identified by this 
method. (See ref. 44.) 


5. Concluding Remarks 


In this paper, we have given the rudiments of generalized function theory and some 
applications in aerodynamics and aeroacoustics. These applications depend on the concept, 
of generalized differentiation and on the Green’s function approach. We have briefly discussed 
the generalized Fourier transformation. Many more examples could be given. The power of 
this theory stems from its operational properties. In addition to the exchange of limit, processes 
that leads to many useful results, discontinuous solutions of linear equations using the Green’s 
function are easily obtained by posing the problem in generalized function space. As seen in thf. 
example of the Oswatitsch integral equation of transonic flow, a nonlinear partial differential 
equation with a discontinuous solution can be cast, into an integral equation based on the 
fundamental solution of the linear part, of the differential equation. The Schwartz generalized 
function theory has unified many ad hoc methods in mathematics and has answered some 
fundamental questions about, linear partial differential equations. The nonlinear theory now 
being developed, in which multiplication of generalized functions is allowed, can be even more 
useful in applications. Generalized function theory is an extension of classical analysis and gives 
engineers and scientists added power in applications. This extension is much like the complex 
analysis that extends real analysis and is very important, in applied mathematics. Finally, 
multidimensional generalized functions, particularly the delta function and its derivatives, are 
quite useful in many applications. 
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